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We propose an inhomogeneous open spin ladder, related to the Kitaev honeycomb model, which 
can be tuned between topological and nontopological phases. In extension of Lieb's theorem, we 
show numerically that the ground state of the spin ladder is either vortex free or vortex full. We 
study the robustness of Majorana end states (MES) which emerge at the boundary between sections 
in different topological phases and show that while the MES in the homogeneous ladder are destroyed 
by single-body perturbations, in the presence of inhomogeneities at least two-body perturbations 
are required to destabilize MES. Furthermore, we prove that x, ^, or z inhomogeneous magnetic 
fields are not able to destroy the topological degeneracy. Finally, we present a trij unction setup 
where MES can be braided. A network of such spin ladders provides thus a promising platform for 
realization and manipulation of MES. 

PACS numbers: 03.67.Lx,71.10.Pm,05.30.Pr,75.10. Jm 



I. INTRODUCTION 

The study of Majorana fermions in various solid-state 
systems has recently attracted a lot of attention-f^Ki^ In 
particular, the possibility of realizing them as zero-energy 
states localized at the end of one-dimensional systems [so- 
called Majorana end states (MES)] has been the subject 
of many recent investigations.^ Besides being of funda- 
mental interest, the study of MES is motivated by their 
potential use for topological quantum computing. 

While electronic systems are a natural choice for 
the realization of such MES, well known fermionization 
techniques have stimulated the study of MES appear- 
ing in the Fermionic mapping of one-dimensional spin 
systems. ^^ ^^ Since many proposals exist for implement- 
ing designer spin-spin interactions,^ and, in particu- 
lar, for the experimental realization of Kitaev-like spin 
models,'^IH^ such one-dimensional systems are interest- 
ing candidates for the realization and detection of MES. 

Proposals of highly engineered spin interactions also 
suggest the implementations of inhomogenous systems 
where regions in different topological phases coexist, 
which is the main feature of the spin ladder studied in 
this work. The spin ladder is based on nearest-neighbor 
Ising interactions, in extension of the Kitaev honeycomb 
model,!^ and consists of sections which alternate between 
the topological and nontopological phases. The phase of 
each section is simply determined by the strength of the 
Ising couplings and at the boundary between different 
phases (as well as at the open ends of the ladder) single 
well- localized MES exist. Being interested in the ground- 
space properties of the model, we perform an extensive 
numerical analysis of possible vortex configurations. We 
always find that the ground state is either vortex-free or 
vortex-full and thereby give strong evidence that that 
the Lieb theorem, originally formulated for different 
boundary conditions, also applies to such open spin lad- 
ders. While freely moving vortices would introduce ad- 
ditional degeneracies, the presence of a finite gap for the 



vortex excitations makes it possible to characterize the 
ground space in terms of spatially localized Majorana 
states, and to consider braiding of such MES. We dis- 
cuss a trijunction setup similar to the one proposed in 
Ref. |15, which could allow the implementation of topo- 
logically protected quantum gates. 

In contrast to previous studies,!^^^!^ we do not analyze 
our model by making use of the Jordan- Wigner transfor- 
mation, but consider the alternative mapping proposed 
by Kitaev.!^ This method has the interesting feature of 
mapping local spin operators to local Fermionic opera- 
tors and has some advantage in analyzing the stability 
of the MES. In particular, it makes immediately clear 
that MES in the homogeneous ladder are already frag- 
ile against single-body perturbations. This conclusion 
is consistent with general arguments showing that topo- 
logical protection in one-dimensional systems can only 
be effective against local perturbations conserving cer- 
tain symmetries.l^^M^In p-wave wires, for example, MES 
are protected against local perturbations which conserve 
parity. Such symmetry assumption is realistic for super- 
conducting wires (although other limitations exist due 
to quasiparticle poisoning^^) but an analogous constraint 
does not appear for general spin systems. For this reason 
we study in great detail the robustness of MES and char- 
acterize the form of local perturbations which can split 
the corresponding degeneracy. In the Fermionic language 
adopted here, the fragility of the topological degeneracy 
results from the only apparent locality of the Kitaev map- 
ping. Non-local strings of Fermionic operators coupling 
the MES app ear through the projector onto the physical 
subspace.^^^ Nevertheless, we show that the inhomoge- 
nous ladder has better properties than the homogenous 
laddei^^ since single-body perturbations are not sufficient 
to split the degeneracy but at least two-body interactions 
are necessary. Furthermore, we show that the topological 
degeneracy of our model is robust against inhomogeneous 
magnetic fields along x, or z. 

The paper is organized as follows. In Sec. [lljwe intro- 
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duce the inhomogeneous ladder and the spin-to-fermion 
mapping we use to study MES properties. In Sec. |III[ 
we derive the topological invariant of the spin ladder for 
both the vortex- free and the vortex-full sectors. We show 
that well-localized MES appear at the junction between 
topological and nontopological sections and demonstrate 
how to move MES along the ladder. In Sec. |IV| we in- 
vestigate numerically the validity of the Lieb theorem^^ 
and show that it also applies to this type of open spin 
ladders, since the ground state is either vortex-free (when 
one of the coupling is negative) or vortex-full (when all 
couplings are positive). In Sec. |v[ we perform an anal- 
ysis of the robustness of MES under local perturbations. 
We demonstrate that single-body perturbations split the 
topological degeneracy of the homogeneous ladder, while 
two-body terms are required in the inhomogeneous case. 
Furthermore, we prove that x, or z inhomogeneous 
magnetic fields alone are not able to split the topological 
degeneracy of the ladder. In Sec. VI we present a trij unc- 
tion setup which could allow braiding of MES following 
the scheme of Ref. [15] and Sec. IVIII contains our final 
remarks. We present additional details and interesting 
aspects of the model in Appendixes A][E 



II. 



INHOMOGENEOUS SPIN LADDER AND 
KITAEV'S MAPPING 
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FIG. 1. Inhomogeneous spin ladder. Each site contains 
a quantum spin-| which interacts with its nearest-neighbor 
spins via bond-dependent Ising interactions Jx, Jy, Jz- The x, 
and z bonds are indicated by red, green, and blue lines, re- 
spectively, the A (B) sublattice site by black (white) dots, and 
the nth unit cell composed of four spins by the black dashed 
square. In contrast to the standard honeycomb model, the z 
couplings are inhomogeneous, i.e., site-dependent Jz Jzij- 



The spin model we propose is an inhomogeneous ladder 
version of the compass^^ or Kitaev honeycomb^ model 
with Hamiltonian 
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where cr^ = (<jf , <j]^, ) are the Pauli operators for the 
spin-1/2 located at the site i of the ladder, and where the 
sum runs over all pairs of nearest-neighbor sites of the 
open ladder containing N unit cells. We assume that the 
ladder is of length 2N—1 (lattice constant set to one) , i.e. 
consists of an odd number 2N — 1 of square plaquettes, 
see Fig. [l] The anisotropy direction in spin space of the 
Ising interaction Jq,.^. depends on the orbital location of 
the bond which is labeled by the index aij = x, 2: for x. 



and z-bonds, respectively, see Fig.[T] Furthermore, we 
allow the z couplings to depend on position, i.e., Jz 
Jz-j. Without loss of generality we assume that Jz^^ > 0. 
Following Ref. IH this model can be solved exactly in 
an extended Hilbert space C by assigning four Major ana 
fermion operators 6^'^'^ and q (all self-adjoint) to each 
site of the lattice and mapping each spin operator to a 
product of two Major anas. 



ibfc, 



(2) 



In fermionic representation, the spin Hamiltonian in 
Eq. (jlj) takes the form H = i AijCiCj^ where Aij = 



ij 



and Uij 



^^cx^j ^(x^j rpj^^ ^ operators commute 
with each other and with iif, and satisfy ufj = 1. There- 
fore, the extended Hilbert space splits into subspaces Cu, 
i.e., C = ^Cui where u represents a certain configuration 
of eigenvalues Uij = ±1. To remove the ambiguity due 
to Uij = —Uji, we assume that for a chosen value Uij the 
first index i belongs to sublattice A (see Fig. [T]). The 
physical subspace is defined through the gauge operators 
Di = bfbybla as £ = {|^) : = 1^)}- Starting 

from a state |^) in the extended space, the corresponding 
physical state |^)phys is given by symmetrization over all 
gauge operators Di, i.e.,^ 



4N 



i^)phys = n 



2 = 1 



1 + A 



(3) 



In each subspace >C^, the operators Aij are replaced by 
numbers Afj and thus the quadratic Hamiltonian H is 
easily solvable with a canonical transformation to new 
Major ana modes^ 



(4) 



Under this transformation, H takes, for a given config- 
uration u, the form Hu = f Em=i ^rn^2m-i^2m, whcrc 
Cm > are the positive eigenvalues of 2iA'^. By defining 
new complex fermion operators = (62m- 1 + ^^2m)/2, 
we finally obtain = ^^(ajn^^m - 1/2). 

The spin ladder Eq. (IT]) possesses as conserved quanti- 
ties two types of plaquettes that are naturally associated 



with each unit cell n, i.e., Wn = 
Wn-- 



"^n,l^n,2<^n,3^n,4 



crJI^i qCtJI /I , where each spin (Jn 



^n,l^n+l,2^n+l,3^n,4 



and 
is la- 
beled by the index n for the unit cell and a = 1, 4 for 
one of the four spins inside the unit cell (see Fig.jl]). We 
say that a square plaquette p carries a vortex if Wn = — 1 
for p = 2n — 1 and if Wn = — 1 for p = 2n. In fermionic 
representation, the plaquettes are products of u opera- 
tors around each unit cell. In the following, we study the 
robustness of MES against local perturbations and show 
that not all the degeneracies due to MES can be split by 
single-body perturbations. Additionally, we show that 
the topological degeneracy cannot be fully split by inho- 
mogeneous X, y, or z magnetic fields. We note that in 
the special case with J^. . = Jz our model is equivalent 
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to the usual honeycomb model studied in Refs. |T4| and 
[3TJ As we show below, the Majorana degeneracy present 
in the homogeneous ladder is less protected in the sense 
that it is fully split by single-body local perturbations. 
Similar homogeneous Kitaev ladders have been studied 
in Refs. |3l] and |32l 



III. TOPOLOGICAL PHASES OF SPIN 
LADDERS 



7i '^y 72 73 74 
1 ^ I \ ^ w ^ \ 



Si 



So 



FIG. 2. Inhomogeneous spin ladder with different topological 
sections. Shown are two topological sections Si and S^ (thin 
z-bonds and J23) separated by a nontopological section 

52 (thick z-bonds Jz2)- For the corresponding Jx,y,z values, 
see main text. The wave functions of the four MES 71,..., 4 are 
mainly localized at the phase boundaries and, for Jx > Jy, on 
the lower sites of the ladder as indicated by the large dots. 

We consider now a spin ladder with sections in differ- 
ent topological phases, 5i, ^2, and ^3, which are distin- 
guished by the value of the Jz-- couplings (see Fig. [2|. 
If we focus on the vortex- free or vortex- full sector, then 
we choose the J^.^. couplings as follows: for the Si and 

53 parts, Jz^^ = Jz^ = Jz3 and IJ^J < \Jx ± Jy\, while 
for the ^2 part we have = Jz2 and 1 1 > \Jx ^ Jy\ 
(see Fig. |2|. Below we derive the conditions for the exis- 
tence of zero-energy MES in the vortex-free (vortex-full) 
sector with the use of the mapping ([2| and prove that 
sections Si and ^3 are topological, while section ^2 is 
nontopological. 

The vortex- free sector corresponds to all Uij = +1. In 
contrast, the configuration where all the u^s along only 
one of the axes (say x axis) take on the value —1 is vortex- 
full. From the explicit expression of A^j we obtain the 
bulk spectrum in the vortex- free sector for Jz- - = J^, 



larger or smaller than one. Therefore, we define the topo- 
logical invariants Ua = — sgn((l — l^f |)(1 — |r2'|)), where 
rf 2 are the two eigenvalues of 7^. The system is in the 
nontopological phase when = +1, and in the topological 
phase with MES when u = —1. From the above explicit 
expression for 7i,2, we obtain the following result for the 
topological invariants in agreement with Refs. il4iand[3T1 



Ui U2 



:sgn(2|J,|-2|J, + J^l). 



(6) 



In the vortex-full sector, the topological invariant is given 
by Eq. (|6| with + Jy replaced with — Jy. From 
Eq. ^ it follows that sections Si and ^3 of our model 
are topological, while section ^2 lies in the nontopologi- 
cal phase. This system thus contains four c Major anas: 
71,4 localized at each end of the ladder and 72,3 at each 
junction between topological and nontopological sections 
of the ladder, see Fig.|2j Note that for > Jy {Jy < Jx)-, 
7i and 73 will reside on the A {B) sublattice, while 72 
and 74 on the B (A) sublattice. For the rest of this 
work we always consider the case where \Jx\ > \Jy\ and 
Jz > 0, since the case \Jx\ < \Jy\ can be treated anal- 
ogously. It is worth pointing out that Ti = T2 when 
\Jx\ = \Jy\^ and consequently all the |(/>n,a=i,...,4| will de- 
crease (increase) with n if both eigenvalues of 7i,2 have 
their absolute values smaller (larger) than one. This ex- 
cludes the presence of MES localized at the right (left) 
end of the ladder. Consequently, Eq. (|6| is strictly valid 
only for \ J^\ ^ \ Jy\. 

The four MES operators 71,..., 4 can be easily expressed 
in terms of Majorana operators q through relation Q: 
71,. ..,4 = T.jQi,,...,4j^j^ where the coefficients Qi^^^^j 
and Qi3(4),j are the elements of the imaginary (real) part 
of the e = eigenvectors of matrix iA'^. Figure [3] shows a 
plot of the wave functions of four Majoranas 71,..., 4 for a 
ladder with N = 60 with all Uij = +1, and Jx = 1.0, 
Jy = —0.4, J^^ = = 0.2 for the Si and ^3 sec- 
tions, and Jz2 = 3 for the middle section 52- We can 
see that the MES 71 and 74 are respectively localized at 
the left and right ends of the ladder, while MES 72 and 
73 are localized at the junctions between topological and 
nontopological sections of the ladder (the precise shape 
of MES can be understood more intuitively by mapping 
the spin ladder to two coupled p-wave superconducting 
wires; see Appendix IA|. The four zero-energy eigenval- 



HX"^^^^ ~ V ^ "^JxJy cos(/c) — (1 — rri)jk^ (5) ues of iAfj (corresponding to 71,. ..,4) reside inside a gap of 



where J 



P 



Ik 



Y^(2 + 2 co^{k)){Jx + Jy^J^^ k is the wave vector, 
and m = 0, 2. 

In the vortex-free sector, the MES eigenvectors (p 
with eigenvalues e = can be shown to satisfy 
the following transfer equations (^^+1,0;+^, 0n,Q;+^)^ = 
Ta{4>n,a-\-^^ 4^n-i,a-\-^)^ (the two labels of (j)n,a corrcspoud 
to the unit cell n and one site a of the unit cell), with 

^ = 0, 2, and Ti 2 = J-^y ^l.y • MES can ex- 

ist only when both eigenvalues of To. have absolute value 



about 1.7 for this choice of parameters . From Eq. (|6| one 
concludes that it is possible to move from the topologi- 
cal to the nontopological phase by changing the relative 
strengths of Jx.y.z- Since MES will exist at the junction 
between sections in different topological phases, MES can 
be created, destroyed, and transported by locally (and 
adiabatically) changing the relative strengths of Jx,y,z 
along the spin ladder. Finally, it is well-known that, in 
principle, exchange interactions can be controlled electri- 
cally (for atomistic or nano structures see, e.g., Refs. [33l 
and [34|) . Thus, applying gates over portions of the spin 
ladder will allow one to move the MES along the ladder. 
This is similar to what is done in superconducting wires 
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FIG. 3. Inhomogeneous spin ladder as defined in Fig. a) 
MES wave functions (corresponding to 71,..., 4) as func- 
tions of site i. The curves for 72,3,4 are shifted vertically for 
clarity. The order used for the site labeling of the spin lad- 
der is shown in b). The circles represent the wave function 
weight of 71 (proportional to the area enclosed by the circle) 
at the corresponding site. For both plots we have all Uij — +1 



(vortex-free), = 60, = 1.0, 



in <Si,3, and = 3 in <S2. 
n = 41 and ends at n = 79. 



-0.4, J,, = J,^ 



0.2 



Section S2 starts at unit cell 



with local tuning of the chemical potential, see Ref. [151 

For the sake of completeness, we study in Appendix [E] 
long-range static correlations which exist in the topolog- 
ical phase only. In the standard honeycomb model spin 
spin correlations vanish rapidly with distance. 



IV. VORTEX-FREE (FULL) GROUND STATE 

In this section we give numerical evidence that the 
ground state is either vortex-free or vortex-full for a cer- 
tain range of Jx.y.z parameters. It is tempting to in- 
voke Lieb's theorenP^ (see also Ref. |37j). However, this 
theorem is not directly applicable to our system since 
it requires periodic boundary conditions in the horizon- 
tal direction or \Jx\ = \Jy\ with open boundaries (when 
the reflection plane is taken to be horizontal and going 
through the middle of the ladder). However, different nu- 
merical checks lead us to conclude that the ground state 
of our spin model is vortex- free/ vortex- full for sgn{Jx) = 
(— /+)sgn(J^) and general J^.^. > 0. Figure H is a plot 
of the single- vortex energy for different N ana couplings 



1.0, J,, 



-0.4, and Jz 



Jz 



J^3 = 0.2 in 



sections Si and tSa, while J^-j = Jz^ = ^ in section 
The vortex energy converges quickly with N and is pos- 
itive. We also see that a vortex lying in the nontopolog- 
ical section ^2 has a larger energy since Jz is stronger 
there. Figure |4] reveals that the energy becomes inde- 
pendent of the vortex position in the bulk of the ladder 
when N is large. However, the ground state is vortex- 
free or vortex-full and such additional degeneracies are 



not present. All ground- state degeneracies can thus be 
associated to MES. 

We have numerically investigated the energy of mul- 
tivortex configurations and found that, although the 
vortex- vortex interaction is attractive, the attraction is 
not strong enough to favor the creation of additional vor- 
tices and the ground state remains a vortex-free state 
(see Appendix [B| for more details). Additional numer- 
ical checks with different Jx,y,z configurations are also 
reported there. For all the numerical checks we per- 
formed, the conclusion remains the same: The ground 
state is vortex- free. Since changing the sign of Jx^y is 
equivalent to u^-^ = — 1 for the corresponding bond, the 
system with ^gii{Jx) = sgn(J^) has thus a vortex- full 
ground state, as expected. Although an analytical proof, 
to the best of our knowledge, is lacking, we conjecture 
that Lieb's theorem can be formally extended to the spin 
ladder considered in this work. 
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FIG. 4. Excitation energy AEy of a single vortex as a func- 
tion of its position p (square plaquette) on the ladder for 
different N. The values of Jx,y and J^.^. are given in the main 
text. The junction between sections Si and S2 is at p = 2N/3 
and between <S2 and S3 at p = 4N/3. Note the slight increase 
of AEv (as compared to the bulk) for N — 45 and 120 when 
the vortex is located at the end of the ladder. 



V. ROBUSTNESS OF THE TOPOLOGICAL 
DEGENERACY 

As mentioned in the Introduction, topological order 
does not exist in one-dimensional systems .^^^^ Espe- 
cially, if no symmetry constraints are present, it is al- 
ways possible to find local perturbations which split the 
Major ana degeneracy. 

In line with this general result, we are thus interested 
in determining what type of local spin perturbations can 
split the topological degeneracy and against which type of 
perturbations the MES of the inhomogeneous ladder are 
robust. Indeed, as we shall see, only certain forms of local 
perturbations can split the topological degeneracy. A 
symmetry-based analysis of perturbations that can split 
the topological degeneracy into different spin ladders has 
recently been performed in Ref. ,38i 
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Here we demonstrate that single-body terms cannot 
fully split the topological degeneracy present in the in- 
homogeneous ladder and that at least two-body pertur- 
bations are required. Furthermore, we also demonstrate 
that a twofold degeneracy always remains when x^y^oi z 
inhomogeneous magnetic fields are applied. The presence 
of inhomogeneities in the ladder is important, since much 
simpler perturbations destroy the topological degeneracy 
of the homogeneous ladder as we show below. 

A. Homogeneous ladder 

Here we first focus on the homogeneous ladder, while 
we study the inhomogeneous ladder in the next section. 
In order to study the robustness of MES, we recall that 
the spin system carries four additional h Majoranas aris- 
ing from the h operators at the ends of the ladder (6f , Ij^^ 
^4Ar-i' ^%n) which are completely decoupled from the 
Hamiltonian Hu. These h Majoranas are always present 
in the model (in the extended space) independent of the 
strength of the couplings Jx,y,z- Such a model carries six 
Majoranas in the topological phase, namely two spatially 
separated c Majoranas and two b Majoranas at each end. 
From Eq. ([2| it is thus clear that a local single-body 
perturbation V = hiaf -\- /i4ArcrJ^ will combine two of 
the three Majoranas at each end of the ladder, leaving 
only one Major ana at the right and one Major ana at the 
left end of the ladder (and thus one zero-energy fermion 
state). It has recently been showiP^that, given a certain 
vortex configuration and lattice topology, only half of 
the states in the extended space are physical and these 
physical states have a definite fermionic parity. If the 
physical states have even parity, then the remaining zero- 
energy fermion state is unfilled, while it is filled for phys- 
ical states with odd parity. In any case this means that 
one of the states of the remaining zero-energy fermion is 
unphysical. We thus conclude that the topological de- 
generacy of the homogeneous ladder has been lifted by x 
single-body terms only. 

B. Inhomogeneous ladder 

Let us now consider the inhomogeneous ladder. The 
ladder possesses eight MES, namely four spatially sep- 
arated 7i (see Fig. [2| and two b Majorana operators at 
each end. An argument identical to the one presented 
for the homogeneous ladder is thus applicable here and 
single-body local perturbations V = hicrf -\- h/^NC^^N ^^^^ 
mix one b with one c Majorana at each end of the lad- 
der. Since four Majoranas are left in the extended space, 
we still have a twofold degeneracy in the physical space. 
The corresponding states in a given u configuration are 
denoted by with m = 1,2. The four Majoranas 

are well- separated from each other; hence, it seems im- 
possible to remove this remaining degeneracy with local 
perturbations V] indeed (^]^|V|^2) = 0- However, one 



has to be careful with the fact that the physical states 
are given by symmetrization over all gauge operators, 
see Eq. ([3|. The physical vortex-free degenerate ground 
states l^if) given by 

\^xi) = \{{}^)m,2), (7) 



where u is the configuration with all Uij = +1. 

From Eq. ([7|), we realize that the projection onto the 
physical subspace introduces a string of Di operators and 
that the spin to fermion mapping ([2| is only apparently 
local. Similar to the case encountered with a Jordan- 
Wigner transformation,!^ the nonlocality of the mapping 
means that a local spin perturbation is nonlocal in the 
fermionic language and can thus connect well-separated 
MES. In order to check whether a local perturbation V 
splits the Majorana degeneracy, we thus need to calculate 
the following matrix elements: 

(vE-n^i^r) = «), (8) 

with m,n G {1,2}. Since V contains string of Di oper- 
ators, we need to check whether V introduces transition 
between 2) and IliGA^*!^! 2)5 where A is a subset of 
indices. If this is the case, then the perturbation V can 
split the remaining degeneracy. 

In the following we explicitly construct a local pertur- 
bation which splits the remaining degeneracy in agree- 
ment with the general result of Refs. l25l - [27l stating that 
topological order is impossible in one-dimensional sys- 
tems. Furthermore, we show that such a perturbation 
has to contain at least two-body terms. 

A necessary criterion for a local perturbation V to split 
the degeneracy is that it does not create vortices. This is 
so because the degenerate ground states are vortex-free 
and states with different u configurations are orthogonal. 
Let us thus first consider all possible single-body terms 
which do not create any vortex and show that they can- 
not split the degeneracy. The effect of Vi = ciaf and 
V4N = ^ANcr^N been considered already. The re- 
maining single-body terms which do not create vortices 
are 

V2 = e2Crf , V4Ar-l = e4Ar-icr|^_i. (9) 

For these single-body perturbations, the only string of 
Di operators which leaves the li-configuration invariant 
is the product of all Di, i.e., Yl^^iDi. As an explicit 
example, let us now consider the effect of V2: 

AN AN 

4llD, = ibyc2l[b^b^btc, (10) 

i i=l 

a iblc2b-,blblj^_,blj^ H U,, H C, (11) 

(ij) ^=1 

(X ib2C27t, (12) 
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where tt = ^2^47v-i^4Ar Ylt^i ^ parity operator and 
we used the fact that ah Uij = +1. We thus have 

{^"Jibyc27r\K) = 7r{<i>-Jiblc2\K) = 0, (13) 

where m,n G {1,2}, tt is the parity of 1^12) and the 
last equahty comes from the fact that C2 creates a finite- 
energy fermion (the contribution of C2 in the remaining 
MES is exponentiahy supressed with system's size). We 
thus finally conclude that single-body perturbations can- 
not fully split the topological degeneracy. 

In the following we show that two-body interactions 
are enough to destroy the remaining MES and we con- 
struct explicitly a perturbation which splits them: 



V 



(14) 



where M is even and belongs to the nontopological sec- 
tion ^2, see Fig. ([2|. An important point to notice is 
that such a perturbation will not create any vortices 
since crM_i changes the values of u\^_2 while 
changes the value of u\^_^ ^. If this was not the case, 
then all matrix elements (^^1^ FIzga ^*l^n) would triv- 
ially vanish, since states with different u configuration 
are orthogonal. The fact that states with different u con- 
figurations are orthogonal implies that the matrix ele- 
ments (^ml riiGA ^^^M-i^Afl^n) can be different than 
zero only if there exists a string of Di operators which 
connect the two different configurations of u. Having in 
mind that Di changes the value of u for the three links 
connected to site z, it is easy to check that Hi^^^ i^" 
troduces transitions between these two configurations of 
u. We have 



M-2 



M-2 



(15) 



2=1 



where the other h operators combine into u operators 
which are all equal to 1. From Eq. (15) we then obtain 



M-2 



M 



^M-l^M 



(16) 



The string of q operators in Eq. ( 16 ) is nonlocal and con- 



nect MES 62 with MES cm- This makes it possible to 



conclude that a?^ icr 



y 

M-l^M 



n 



M-2 



Di can split the remain- 



ing topological degeneracy. However, checking that the 
matrix elements are really different from zero, i.e.. 



M-2 



(17) 



requires determining whether the string operator fli^i 
makes it possible to come back to the ground subspace. 
This is, in principle, the case since q's are superposi- 
tions of all eigenmodes. In Appendix [C] we use a different 
approach and show explicilty that V can indeed induce 
splitting between the Majorana states (all these consid- 
erations remain valid for a perturbation V = ea^a^^^ 



with odd M as well). Although single-body terms are not 
able to fully split the topological degeneracy, we showed 
that local two-body terms are. From this result, it ap- 
pears that MES are quite fragile against external mag- 
netic fields since two-body interactions are easily gener- 
ated in second-order perturbation theory. However, on 
the positive side, we show below that the topological de- 
generacy is protected against inhomogeneous x, and z 
magnetic field components. 



C. Protection against magnetic fields aligned along 
X, or z direction 

Let us consider a local perturbation V which represents 
an inhomogeneous magnetic field in direction a = x^y^z: 



V = J2hiar. 



(18) 



As we demonstrated in the previous section, single-body 
terms cannot fully split the topological degeneracy. As a 
consequence, the perturbation V is not able to split the 
topological degeneracy at first order. We thus ask here 
the question whether this is possible at higher orders. 
The answer, in principle, is yes. For MES well-separated 
by L sites, the degeneracy can obviously be split in L*^ 
order. However, L*^-order terms are exponentially sup- 
pressed with system's size and the effect of such pertur- 
bation is thus negligible for large systems. We are then 
interested in determining whether this is possible at order 
n independent of system's size. 

As mentioned previously, a necessary condition for a lo- 
cal perturbation to split the degeneracy is that it does not 
create any vortices. The good candidates V"^ generated 
at order n have thus to be products of terms and 
^1 AN ^2 4 AT- 1 5 where z,j are nearest-neighbor sites 
along an a-link (i, j)^. Note that of this form not 
only do not create vortices, but also do not change the 
configuration of u. Thus, the two matrix elements which 
need to be considered are 

{^-JV-\K) and {^-JV-l[Di\^l), (19) 



with m, n G {1,2}, since Y[t=i leaves the configuration 
of u invariant. If we recall that (the remaining four) MES 
are well- separated from each other, it is then clear that 
= for n < L. Furthermore, as mentioned 
previously, the operator fli^i A is proportional to the 
parity operator of the eigenmodes^^ and we thus have for 
n < L 



4N 



i^-jv" n ^^1*") ^ ^C^liv-IK) = 0, (20) 

i=l 

where tt is the eigenmodes parity of |^i^2)- 

From this we conclude that the effect of x, or z 
magnetic field components on the inhomogeneous spin 
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ladder is exponentially suppressed with system size. In 
other words, the topological degeneracy cannot be fully 
split by inhomogeneous magnetic fields purely along the 
X, or z direction. 



VI. BRAIDING MES IN A TRIJUNCTION 
SETUP 

The recent demonstration of non-Abelian character of 
MES in p-wave wires^^ makes the possibility to realize 
topological quantum computing by braiding MES an at- 
tractive and promising method. Majorana end states in 
spin ladders can be moved by locally tuning the value 
of Jz couplings (similar to the local tuning of chemical 
potential in p-wave wires) and it is thus in principle also 
possible to braid them (following the scheme of Ref. |15|) 
in the trijunction setup presented in Fig. [5] The trijunc- 
tion in Fig. [5] possesses the interesting property that no 
spurious Majorana modes are created in the course of 
braiding. Indeed, when all three parts building the tri- 
junction are in the topological phase, the Jx^y couplings 
(dashed lines in Fig. 5) between the MES will combine 
two of them into an ordinary complex fermion. We note 
that such braiding processes in the simpler xx — yy chain 
would be difficult to realize because of the many addi- 
tional zero-energy modes, see Appendix [D] This is why 
we are interested in inhomogeneous spin ladders where 
such additional degeneracies do not appear and all Ma- 
jorana modes are well-separated from each other. 

In order to be of use for topological quantum comput- 
ing, an important property of MES is that they follow 
non-Abelian exchange statistics (as it is the case in p- 
wave wireJ^. As mentioned above, the projection op- 
erator (onto the physical subspace) V introduces strings 
of operators and renders mapping ^ only apparently 
local. In such a case the study of braiding statistics is 
more complicated and especially mapping the Hamilto- 
nian onto Kitaev's toy model Hamiltonian is not enough 
to conclude anything about the statistics of MES.^^ We 
leave here the question of the MES statistics as an open 
problem. 

To perform a robust topological quantum computa- 
tion by exchanging MES, the two-body perturbations 
splitting the topological degeneracy have to be excluded. 
This requirement is analogous to the parity conserva- 
tion requirement in p-wave wires and can, in principle, 
be achieved by carefully screening magnetic fields. Since 
Isin g inte ractions can be tuned by means of electric fields 
QY^ly^ESSM ]2iagnetic fields are not required for braiding 
and there is no intrinsic contradiction between exchang- 
ing MES and screening external magnetic fields. 



VII. CONCLUSIONS 

We have proposed inhomogeneous spin ladders and 
shown that they support a topological phase with local- 




FIG. 5. Network of spin ladders, a) Majoranas are exchanged 
through the trijunctions, which are shown in detail in b). The 
connection between the three spin ladders is given by Ising 
couplings Jx (red dashed lines) and Jy (green dashed lines). 
Braiding is performed by varying J^, see main text. MES 
7i,...,4 (large dots) are localized at the left and right ends 
of the ladder and at the junction between topological (thin 
2;-links) and nontopological (thick z-links) sections. 

ized Majorana states. We have studied the robustness of 
MES under local perturbations and demonstrated that 
single-body perturbations are not enough to split the 
topological degeneracy. We have explicitly constructed 
a two-body perturbation which does this, in agreement 
with the general fact that topological order is not possi- 
ble in one-dimensional systems.^^^^^On the positive side, 
we showed that the topological degeneracy cannot be de- 
stroyed by inhomogeneous magnetic fields aligned purely 
along X, y^ or z direction. Finally, we have presented a 
trijunction setup where MES can be exchanged similar 
to Ref. 15, 

While the spin ladders envisioned here are not yet 
available, we hope that the present study provides a 
strong encouragement for their experimental realization 
since they represent a promising framework for the real- 
ization and manipulation of MES. 
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Appendix A: Mapping to two coupled Kitaev p-wave 
superconducting wires 

As presented in the main text, the model we consider 
possesses three different sections 5i, ^2, and ^3 which 
are distinguished by the value of the Jz^j couplings [see 
Fig. [g] . We focus on the vortex-free and vortex-full sec- 
tors where we choose Jx,y,z,z' such that Si and ^3 are 
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71 



72 
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74 



^1 



FIG. 6. Inhomogeneous Kitaev spin ladder. This spin ladder 
possesses two topological sections Si and S3 (thin z links with 
couplings Jzi and Jz^ which we choose to be equal, i.e., Jzi = 
J23 — Jz) separated by a nontopological section <S2 (thick z 



links with couplings J^. 



Jz')- The main components of 



the four MES wavefunctions 71,..., 4 lie on the lower sites for 
Jx > Jy and are represented by large dots. 



topological, while section S2 is nontopological. This sys- 
tem carries four MES: 71 and 74 at the left and right end 
of the ladder, respectively, while 72 and 73 sit at the junc- 
tion between topological (52) and nontopological (Si^s) 
sections of the ladder. 

Let us focus on a topological section, say 5i, in the 
vortex- free sector (i.e., Uij = +1), and study the location 
of MES 7i,2 and their behavior under the modification 
of Jx,y,z couplings. We disregard here the presence of 
the two other sections 52,3- It is useful to consider our 
spin system as two xx-yy chains coupled via Jz Ising 
couplings. Let us now introduce the following complex 
fermion operators: 



1 



(c2j-i ^ic2j) and d\ = ^(c2j-i - ic2j), (Al) 



with j = 1...2N (the site labeling is shown in Fig. [3|, 
{dj^dj'} = 0, and {dj^d^j,} = djj'. Then the upper (u) 
xx-yy chain is mapped to the Kit aev mo del for a one- 
dimensional p-wave superconductor 



N N-1 



H"" = -^""^dldj-Y^ (t''d]dj^i^A''djdj^i^h,c}j , 



(A2) 

with /i^ = 2Jx and = — = J^, while the lower (/) 
xx-yy spin chain is mapped to 



2N 



2N-1 



= -jjL^ ^ d]dj- ^ [t^ d]dj j^i + A^djdj^i + h.c. 
j=Ar+i j^N-\-l 

(A3) 

with /i^ = 2Jy and = — = J^. 

The Jz spin couplings between upper and lower xx-yy 
chain leads to a hopping term H^^ between upper and 
lower wire in the fermionic representation, 

N 

H^' = -Y, {t^'d]d2N-U-i) + h.c.) , (A4) 



where t^^ = 2Jz- 

Let us first focus on the case = 0. Then, the system 
consists of two decoupled wires, Eqs. (A2) and (A3), 



and we can distinguish between the following cases: If 
Jx > Jyi then the upper wire lies in the nontopological 
and the lower wire in the topological phase, i.e., the two 
MES are localized in the lower wire, one at the left and 
one at the right end; vice versa for Jy > J^. 

When the z-couplings are turned on, i.e. Jz > 0, then 
the MES spread over both the upper and the lower wires 
as shown in Fig. 3 of the main text. If Jz increases, 
then the MES continue to spread until they completely 
split when \Jz\ > \Jx -\- Jy\ in the vortex-free sector and 
\Jz\ > \Jx — Jy\ in the vortex- full sector, see Eq. (6) in 
the main text. It is also straightforward to understand 
the exact site localization of the two MES. For J^ > Jy, 
most of the weight of the left 71 (right 72) MES resides at 
respectively the first and last site of the lower xx-yy chain 
and spreads only over A {B) sublattice sites. Indeed, 
the Ja;^^ ^2 -couplings between spins residing on different 
sublattices forbids 71 (72) to spread over B (A) sites. 
Similarly, for Jy > J^, most of the weight of the left 71 
(right 72) MES resides at, respectively, the first and last 
sites of the upper xx-yy chain and spreads only over B 
{A) sublattice sites. 



Appendix B: Vortex- free and vortex-full ground 
states 



As discus sed in the main text, although Lieb's 
theorenJS^ is not directly applicable to our system, we 
nevertheless are able to show numerically that the ground 
state is indeed vortex- free for sgn(Ja^) = —sgii{Jy) and 
Jz^j > 0, while it is vortex-full for sgn(Ja^) = sgn(J^). 
Let us focus on the case sgn{Jx) = — sgn(J^), since the 
other one can easily be deduced from it as discussed in 
the main text. Figure [7| shows single- vortex energies for 
different ladder lengths and Jx,y,z,z' coupling configura- 
tions. All the results are consistent with our assumption 
that the ground state is vortex-free. We have further- 
more investigated the effect of vortex- vortex interactions 
and plotted in Figs. [S] and [9] mult i vortex energies for dif- 
ferent TV, Jx,y,z couplings, and distance between the vor- 
tices. Again, all the plots indicate a vortex-free ground 
state since the attractive vortex- vortex interaction is not 
strong enough to favor the creation of additional vortices. 
' Finally, we plot in Fig. [lO] the energy of the vortex-full 
sector as a function of N for J^ = 1.0, Jy = —0.55, 
Jz =0.25 in 5i,3, while J^/ = 4 in 52- The energy of the 
vortex-full sector is always positive and increases linearly 
with the system size N. This result again shows that 
vortex-vortex interactions do not favor the creation of 
vortices and the ground state is free of vortices. We have 
checked that this result is valid for many other choices of 
parameters Jx,y,z- A detailed explanation of the plots is 
given in the figure captions. 

As a final remark, we would like to mention that the 
groundstate of the trij unction setup presented in the 
main text (see Fig. 5 of the main text) is also vortex-free 
or -full. Indeed, each ladder forming the trijunction is 
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separately free or full of vortices and thus by a continuity 
argument it is clear that switching on (small) couplings 
between different ladders cannot create vortices. 
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FIG. 7. Energy AEy of a single vortex as a function of its 
position p on the ladder. We recall that a vortex can be placed 
at 2N — 1 different positions on a ladder with N unit cells. 
The five different curves correspond to N = 9, 15, 30, 60, 150. 
We see a clear difference between the vortex energy in the 
bulk and near the boundaries: Boundary effects increase the 
energy of a vortex lying near to one end of the ladder. It 
is also worth pointing out that the vortex energy converges 
quickly (with N) to its thermodynamic limit value. We see 
that the vortex energy is positive for each curve irrespective of 
the vortex's position. This plot thus supports our claim that 
the ground state is vortex free. The values of the different 
couplings are: a) Jx = 1.0, Jy = —0.5, Jz = 0.3 in <Si,3, and 
J^/ = 0.3 in ^2. h) Jx = 1.0, Jy = -0.65, Jz = 4.3 in 5i,3 
and Jz' = 4.3 in S2 ■ The curves for N = 15, 30, 60, 150 are 
shifted vertically by 0.005,0.01,0.015, and 0.02, respectively, 
for clarity. 



Appendix C: Different mapping to study the 
robustness of MES 

The aim of this Appendix is to propose another map- 
ping to study the robustness of MES in the homogeneous 
and inhomogeneous ladder. We use the same labeling as 
Fig. [3] but now make use of the following spin-to-spin 




0.005 



250 



FIG. 8. a, 6) Energy AEv of a single vortex as a function of 
its position p on a ladder with N = 120, Jx = 1.0, Jy = —0.37, 
Jz = 0.25 in tSi,3, while Jz' = 4 in <S2 for a) and Jx = 1.0, 
Jy = -0.37, Jz = 0.25 in <Si,3 and Jz' = 0.25 in S2 for b). 
The junction between sections <Si and S2 is at p = 2N/3 and 
between S2 and S3 is at p = 4:N/3 



mapping 



mm 



i—1 ' i "> 



4Ar 

n 



k 5 



(Cl) 
(C2) 



where r^^^'^'^ are usual Pauli matrices. The mapping is 
thus to an auxiliary system of 4A^ + 1 spins (i = 
0, 1, 2, 4A/'). Therefore, the mapping introduces a dou- 
bling of the Hilbert space which can be taken into account 
by imposing the constraint Tq = +1. 

The spin ladder Hamiltonian ([T]) takes in this new lan- 
guage the following form: 

R = J^{t^^t1^...^t'^_^ 



+ 1 



I T ( z y X y z 



Ar-4'^Ar-3''A/'-2''Ar 



y 



(C3) 

This Hamiltonian possesses many conserved quantities 
since all with odd i commute with H and are straight- 
forwardly connected to the vortex operators: 



'2n-l '2n+l 



and Wr, 



"'^2n+l'^2n+3- i^^) 
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FIG. 9. a, 6) Energy AE2V of two vortices as a function of 
the position p of the second vortex. The first vortex lies on 
the p = 1 square plaquette. The Jx,y,z,z' parameters are cho- 
sen respectively as in Fig. [S] The vortex- vortex interaction 
is attractive and rapidly decaying as a function of distance 
between the two vortices. Indeed, already for p = 5 the en- 
ergy of the two vortices is roughly 0.04 [0.03 (energy of the 
vortex at the boundary p = 1) plus 0.01 (energy of the vor- 
tex in the bulk)]. However, as mentioned in the main text, 
the attraction is never strong enough to favor the creation of 
vortices and the energy of the two vortices is always positive. 
The junction between sections Si and S2 is at p = 2N/3 and 
between S2 and S3 is at p = 4:N/ 3 



FIG. 10. Energy Evf of the vortex- full sector as a function 
of N with Ja^ = 1.0, Jy = -0.55, Jz = 0.25 in <Si,3, while 
Jzf = 4 in <S2. As expected, the energy of the vortex- full sec- 
tor is always positive and grows linearly with N. The slope of 
the of the straight line can be interpreted as an average vortex 
energy. This plot indicates again that the vortex- vortex inter- 
action does not favor the creation of vortices and the ground 
state is vortex- free. The junction between sections Si and & 
is at p = 2N/3 and between S2 and S3 is at p = 4A/'/3. 



In the strong Jx limit, the two corresponding degenerate 
ground states take the following form: 



l^n) = I to +2 +4 ....+Ar-2 Un), 



(C5) 



where +i is an eigenstate of rf while tj and are eigen- 
states of Tj. 

It is now straightforward to understand that the single- 



body perturbation V 
cal degeneracy since 



— ^^^^ ^P^i^ topologi- 



(V'tl^l^4) = (V'tkivlV'i)/o. 



(C6) 



With the Kitaev's mapping, we noticed the presence of 
an additional degeneracy because of the presence of six 
MES in the extended space. In the new language, this 
additional degeneracy corresponds to the sign of rf (the 
state with all T2^_i = +1 is degenerate with the state 
where all t|^_i = —1). This degeneracy is split by af = 

As a first step, let us now study the homogeneous lad- Iljli since it induces transition between the states 



1. Homogeneous ladder 



der and consider the case where Jx ^ Jy,z which, accord- 
ing to Eq. (|6|, is in the topological phase. For Jy^^ = 0^ 
the odd spins are completely decoupled and the corre- 
sponding 2^/^~^ degeneracy is due to the absence of gap 
to create a vortex. However, there is an additional and 
interesting degeneracy coming from the fact that the last 
spin does not appear in the Hamiltonian and the total 
degeneracy in the topological phase is thus 2^/^. In the 
opposite limit, i.e., when 7^ and Jx^y = 0, which, 
according to Eq. (|6|, corresponds to the nontopological 
phase, then by construction Tq = +1 and no additional 
degeneracy is present. In this limit the degeneracy is 
thus 2^/^~^. We thus identify the additional degener- 
acy due to the N^^ spin as the topological degeneracy. 



with all r. 



2i-l 



+1 and all 



= -1. 



2. Inhomogeneous ladder 

Let us now consider the inhomogeneous ladder with 
two topological sections Si^s separated by the nontopo- 
logical section see Fig. |2l In the case where Jx 
and Jy^z = in 51,3, while Jz ^ and Jx^y = in 52 it 
is easy to see the origin of the eightfold degeneracy, see 
Fig. [11] One degeneracy is related to the sign of rf = ±1. 
This degeneracy is easy to split with a local perturbation 
[as we showed with Kitaev's mapping ([2|, 
f also splits this degeneracy]. The origin 



err 



nN 



11 



where Wx = = Jx/^ and Wy = — = —Jy/A. 

Since there is a difference of tt in the pase of A^^ and 



A^, we thus conclude that Hamiltonian (D3) represents 



of the remaining fourfold degeneracy can be understood 
easily in the limiting case we consider. The first twofold 
degeneracy is due to the N^^ spin which is decoupled 
from the other spins while the other twofold degeneracy 
comes from the middle Ising section, which has the same 
energy if all spins are up or down. The former degener- 
acy is trivially split by a perturbation erfj = ea^, while of freedom and rewrite Eq. (|D3|) as 



an array of tt junctions and thus possesses 2N additional 
zero-energy modes. To find the spectrum of Hamil- 
tonian (D3), we artificially double the number of degrees 



the later is clearly split by the perturbation in Eq. (14), 



(C7) 
(C8) 



where a 



ai 



t 



a2N CLi 



(D4) 



and H is a real 



since r^_^ is a conserved quantity (M - 1 is an odd site 47V x 47V symmetric matrix defined through Eq. ( [D3 ). 
here). We have thus shown here explicitly that two-body 
perturbations are enough to split the remaining degen- 
eracies in agreement with the discussion in the main text. 



Figure [12] is a plot of the eigenvalues oil-L which corre- 
sponds to the modes of a xx-yy spin chain with Jx = 0.4, 
Jy = 1.0, and length 27V = 20 for a) and 27V = 100 for 
b). As expected, the number of additional zero-energy 
modes is 27V. It is thus possible to generate zero-energy 
modes in the xx-yy spin chain by increasing the system 
size. 



FIG. 11. Pictorial representation of the groundstates in the 
limit Jx and Jy^z — in sections <Si,3 and 7^ and 
Jx,y = in 82- Here ft and ^ are a pictorial representations 
of eigenstates. 



Appendix D: Proliferation of vr-junction zero-modes 
in XX-YY spin chain 

In this section we study some properties of the zero- 
energy modes present in a xx-yy spin chain of length 27V, 
described by the Hamiltonian 



XX— yy 



2N-1 

7 \ ^ ^X 

iodd 



2N 



y^y 



(J, (J. 



+1- 



(Dl) 



We show that the xx-yy spin chains contain many ad- 
ditional zero-energy modes besides the two c Major anas 
localized at the ends of the chain. In the language of map- 
ping this arises because all non equivalent u config- 
urations are degenerate; that is, putting a i^^- = — 1 does 
not cost energy. It is instructive to study this model with 
a Jordan- Wigner transformation 



/c=l 



where annihilates a complex fermion at site j, i.e., 
{a^-^\a^t^} = and {aj,a|/} = ^j^/, and 



With the use of Eq. (D2), i^, 



33 

takes the form 



Hxx-yy = ^ (^-Wxalai^i + Axttitti^i + h.c.^ 



i odd 

E 



{-Wyalai+i + Ayttitti+i + h.c.) , 



(D3) 




FIG. 12. Energy eigenvalues Ck of 1-L in Eq. (|D4| for xx-yy 



chains with 0.4, Jy 1.0, and 27V 20 foro) and 27V 
100 for b). There are 27V zero-energy modes in addition to 
the two expected MES. The presence of these additional zero 
modes can be understood by mapping Hamiltonian Hxx-yy 
to an array of tt junctions [see Eq. |D3|. 



Appendix E: Long-distance spin-spin correlation 
function 



In this section we study the static long-distance spin- 
spin correlation function (crfaj^) (the site labeling is 
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shown in Fig. 3 of the main text). We note that this cor- 
relator vanishes in the standard honeycomb modeP^^ 
but is non zero for the ladder in the topological phase 
due to the presence of MES localized at sites 1 and AN 



when Jx > Jy (the scenario with < Jy can be treated 
analogously by considering {o'2crlj^_-^)). Let us first give 
an explicit expression for {crfcr^N)- 



Since 



(El) 



and 



(6i, 62, ^4Ar-i, ^4Ar)Q^ = (ci, C4Ar), 



(E2) 



we have 



N 



(E3) 



Using 



and 



c^CJ=Y,Qk^QlJhbk'' (E4) 

k,k' 



4 = (^2/e-l - ^^2/c)/2 

dk = {hk-i +i62/c)/2, (E5) 



= ^ Q2k-liQ2l-ljb2k-lb2l-l + ^ Q2k-liQ2ljb'2k-lb2l + ^ Q2feiQ2Z-lj^2/c^2/-l + ^ Q2kiQ2ljb2khl , (E6) 



we obtain 



CiCj = Qlk-iiQli-iMk + al){ai + a|) + ^ Qlk-uQliMk + 4)(lA)(a/ - a^^) 

l,k l,k 

+ EQ2fciQ2;-ii(lA)(«fe - + EQ2;ciQ2y(lA')(«fe - - «i) > (E7) 



and thus 



^i^j = [Q2k-iiQ2i-ij{o^kai + aka] + + ala]) 

k,l 

^Q2k-iiQ2iji^/^){^kai - aka] + alai - ala]) 
-^Q2kiQ2i-ij{^/^){(^kai + aka} - a^a^ - ala}) 

^Q2kiQ2ij{-^kai + «/ctt| + alai - 44^) . (E8) 

It is now straightforward to calculate (^n=o|c"i ^''fArl^n^o), where n (n = 0, 1) represents the filling of MES while all 
the high-energy modes are unfilled, 

(^n=o|criCrJ^|^n=o) = -m^4Ar(^ = 0|ciC4Ar|n = 0) 

= -K4Ar(^ = ^\Y1 [Q2fe-liQ2fe-li(^^2fe-l4/c-l + 4/c-1^2fe-l) + Q2fe-liQ2fej (1 A) (-«/c4 + 4^^) 

+Q2feiQ2fc-y (l/0(«fe4 - 4«/^) + Q2fc,Q2fei(«fe4 + aW)] \n = 0), (E9) 
where i = 1 and j = AN. 
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With the use of the fermionic anticommutation relation {a^^al,} = 1 we obtain 

{^n=oWi(^4N\^n=o) = "^^MAT ^ [Q2k-liQ2k-lj " (1 A)Q2/c-liQ2/cj + (1 A)Q2/ciQ2/c-lj + Q2/ciQ2/cj] 

k 

+ J2 [{Vi)Qlk-iiQlk, - (2A)g2fe,Q2fc-y ] {nk = 0) (ElO) 

k 

= -iul,N E [Q2k-uQ2k-i, - (iA)Q2fc-iiQ2fci + (iA)Q2fciQ2fc-ii + QlkiQlk,] ■ (eh) 

k 

Since the matrix is orthogonal we finally obtain 

{^n=oWfaU^„=o) = -iul^N E [-(l/0Q2fc-HQ2fci + (lA)Q2fciQ2fe-y ] (E12) 

= ^lAN ^ [Q2/c-liQ2/cj ~ Q2/ciQ2/c-lj] •> (E13) 

where we recall that u\^^ = ±1 decouples from the Hamiltonian in the absence of external perturbations. 
Similarly we can show that 

k 



where a is the index of the fermonic mode formed by the 
Majoranas; that is, tIq, = n = 1 is the filling of MES. As 
mentioned above, the long-distance spin-spin correlator! 
{^i^In) vanishes in the standard honeycomb modeP^^ 
and is non zero here only in the topological phase due to 
the presence of the MES state with components on both 
sites 1 and 4.N . We show in Figs. [Tsja) and h) a plot of 
— {^nWi^tN\^n) as a function of N with all Uij = +1, 
Jx = 1.0, Jy = -0.4, and J^^ J^^ Jz3 0.2 for 
the topological phase in a) and J^^ = = Jz3 = 2 
for the nontopological phase in b). We make use of the 
projection protocol of Ref.[29lin order to determine if the 
physical ground state of the vortex-free sector has even 
(n = 0) or odd (n = 1) parity. As expected, the long- 
distance spin-spin correlation takes a finite value in a) 
while it vanishes in b). 

In the remainder of this section, we want to investigate 
the effects of vortices on the long-distance correlation 
function. Figure 14 1) shows a plot of — (^n|c^i c^JatI^u) 



as a function of position of a single vortex, p, for A/" = 50, 
= 1.0, Jy = -0.4, J^^ = = Jzs = 0.2. This ladder 
has one topological section with two MES 71,2 localized 
on the left and right ends. The oscillations between pos- 
itive and negative values of the correlator show that the 
vortex changes the value of (^n|c'"i c'"4Ar|^n) and thus the 
MES parity ^71 72 as a function of its po sition on the lad- 
der. Indeed, using Eqs. (E13) and (E14) we show numeri- 
cally that (^n=o|c^fc^J^[l^^o) = -=T^=1 (jj^ |^n=l) , 

and thus conclude that a change of sign in the correlator 
implies a change of the parity ^71 72 (ie.,n = 0^n = l). 
We make use of the projection protocol of Ref. I29> in or- 
der to determine if the physical ground state of the one- 
vortex sectors have even (n = 0) or odd (n = 1) parity. 
In Fig. 14 b) we plot — (^nlcrf crfArl^n) as a function of 




FIG. 13. -(^n|cricrjjv|^n) as a function of A, with all mj = 
+ 1, Jx — 1.0, Jy — —0.4, and Jz^ — Jz2 — Jz3 — 0.2 for a), 
and Jzi = Jz2 — Jz3 = 2 for b). We make use of the projection 
protocol of Ref. 29 in order to determine if the physical ground 
state of the vortex-free sector has even (n = 0) or odd (n = 1) 
parity. 



position of a single vortex, p. 



Jy 



-0.4, J,, = J, 



0.2, J, 



for A" = 50, = 1.0, 
= 2. This ladder carries 



four MES: 71,4 at respectively the left and right ends of 
the ladder and 72,3 at the junction between topological 
and nontopological sections. The oscillations in the cor- 
relator demonstrate again oscillations in the parity ^71 74. 
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;i;ii!i!ii;iiii;:iii!ii;iiii;;:ii:ii;iih;ii:i 



As mentioned in the main text, the one- vortex state is 
highly degenerate since it does not cost energy to move 
a vortex to a nearby plaquette, and thus, without any 
prior measurement, the position of a vortex is generahy 
not known, and so neither is the parity of the MES. In 
the case where MES are used for topological computing 
(braiding), it is thus essential that the ground state is 
vortex-free or vortex-full. 



:iiii'ninnni!i;i!l'I'!!'!I!'I!li!lf!i!'l'' 



FIG. 14. Plot of correlator -{"^nlcrfa^Nl^n) as a function of 
position of a single vortex, p, for N = 50, Jx = 1-0, Jy = —0.4, 
and Jzi — Jz2 — — 0.2 in a) and Jz^ — Jzs — 0.2, — 2 
in b). The junctions between sections <Si,3 and S2 are at pla- 
quettes p = 41, 71. We used the projection protocol of Ref.[29l 
to determine if the physical ground state of the corresponding 
single- vortex sector has even (n = 0) or odd (n = 1) parity. 
Note that the physical groundstates we consider have a fixed 
parity 27273 = +1 and oscillating parity 27174. 
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